Quantum Dynamics of the Avian Compass 
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The ability of migratory birds to orient relative to the Earth's magnetic field is believed to involve 
a coherent superposition of two spin states of a radical electron pair. However, the mechanism 
by which this coherence can be maintained in the face of strong interactions with the cellular 
environment has remained unclear. This Letter addresses the problem of decoherence between two 
electron spins due to hyperfine interaction with a bath of spin 1/2 nuclei. Conditions necessary 
for long lived coherence are identified, and a simple yet robust model for sensing magnetic field 
orientation is presented. 



The ability of a migratory bird to orient itself relative 
to the Earth's magnetic field is at once a familiar feature 
of everyday life and a puzzling problem of quantum me- 
chanics. That birds have this ability is well established 
by a long series of behavioral experiments. However, the 
precise mechanism by which an organism may sense the 
orientation of the weak geomagnetic field remains unclear 
and theoretically problematic. 

Although commonly referred to as the "avian com- 
pass," the ability to sense the local magnetic field ori- 
entation has been observed in every major group of ver- 
tebrates, as well as crustaceans, insects, and a species of 
mollusc 0, 0- For the majority of species, the primary 
compass mechanism appears to be light-activated, with 
only a few exceptions such as the sea turtle or the subter- 
ranean mole rat[2|. In addition to a light-activated com- 
pass located in the eye, migratory birds are believed to 
possess a seperate mechanism involving magnetite, with 
possible receptors identified in the beak[3], the middle 
earQ and the brain stemjB)]. This paper addresses the 
light- activated mechanism, which in addition to being 
widespread is also well studied by a long series of behav- 
ioral experiments, reviewed in 0,0-0] . 

The basic parameters of the compass mechanism can 
be mapped out by confining a migratory bird in a coni- 
cal cage during the preferred migration time. The rest- 
less nocturnal hopping behavior, or Zugunruhe, will then 
tend to orient itself in the preferred migratory direction. 
The importance of various environmental features can 
then be tested by observing their effects on the bird's 
ability to orient. Experiments along these lines have 
shown that the primary compass mechanism is light ac- 
tivated, with monochromatic 424 nm (blue) and 565 nm 
(green) light sufficient to orient, but with an abrupt 
transition to disoriented behavior between 560.5 nm and 
567.5 nm. Under low irradiance 617.0 nm red light, birds 
showed oriented behavior rotated ±90° relative to the 
preferred orientation, behavior which was also observed 
at high irradiance with 565 nm light. [l|, 0] This has 
been interpreted as indicating multiple light-activated re- 
ceptors, with a high efficiency, short wavelength recep- 
tor yielding the correct orientation and a low efficiency, 



long wavelength receptor yielding the rotated orientation. 
Doubling the magnetic field causes temporary disorien- 
tation lasting approximately an hour but correct orienta- 
tion thereafter [lOJ, showing that the compass mechanism 
can adapt to various magnetic field strengths, but that 
this adaptation is not instantaneous. 

Experiments with oscillating magnetic fields have 
shed light upon the quantum mechanical nature of the 
avian compass. Here the birds were subjected to a weak 
oscillatory field oriented transverse to the static field. 
When the frequency of oscillation was narrowly tuned 
to the Larmor frequency for an electron in a static field 
to flip its spin, the birds became disoriented, even when 
the oscillatory field strength (15 nT) was extremely weak 
relative to the static field strength (46 /iT). When the 
static field strength was doubled, so too was the fre- 
quency needed to disorient the birds. These results have 
been interpreted as lending support to a "radical pair" 
model, in which the field orientation is sensed via its ef- 
fect upon two unpaired electron spins. 

As originally formulated in [12, EH, the radical pair 
model holds that the absorption of a photon creates an 
excited state of two unpaired electrons initially in the 
s, m s ) = |0,0) singlet state \s) (or, with identical logic, 
the \s,m s ) ~ |1,0) triplet state \t)). Interaction with 
the local magnetic field then drives coherent oscillation 
between the singlet and triplet states with frequency lo = 
B z Afi. If the singlet and triplet states react at rates k s 
and k t to form distinguishable byproducts, the detection 
of these byproducts constitutes a probe of the quantum 
state [IH and, indirectly, of the field. If k = k 8 + kt, the 
fraction of triplet byproducts is given by 



triplet fraction 



(1) 



While the radical pair model does an excellent job of 
explaining the broad features of behavioral experiments, 
it has both practical and theoretical difficulties. As a 
practical problem, the control signal given by Eq. [T] is 
weak unless k and u are closely balanced. However, dou- 
bling the local field strength yields only temporary dis- 
orientation, while birds in the wild can adapt to local 
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field strengths ranging from 30 /zT at the equator to 60 
/iT at the poles [15 1. 

A more fundamental difficulty with the radical pair 
model is the need to maintain spin coherence in an en- 
vironment which is extremely hostile to it. In a cellular 
environment, every water molecule possesses two spin 1 /2 
protons, which interact with unpaired electrons via the 
hyperfine interaction. However, as observed in [l6j |. co- 
herence between electrons must be maintained for 10~ 4 s 
or longer, or else the rate of spin flips due to an oscil- 
latory field in would be too low to affect the signal. 
Thus, understanding the avian compass requires an un- 
derstanding of decoherence induced by the hyperfine in- 
teraction. In [l7j |. this interaction was converted to an 



of the form 



effective magnetic field, while [14( considered the effects 
of spin relaxation on the radical pair model. [3] propa- 
gated a system of two electrons interacting with several 
nuclear spins, finding that coherence between the elec- 
trons enhanced the sensitivity to magnetic field orienta- 
tion. The related problem of electron spins in a quantum 
dot has been studied in lij 2oj |. 

This paper addresses the problem of a two spin sytem 
in the m s — subspace, where decoherence is induced 
by hyperfine interactions with a bath of spin 1/2 nuclei. 
The form of this decoherence is found to differ if the elec- 
trons interact primarily with the same nuclei (proximate 
electrons) or with different sets of nuclei (distant elec- 
trons), or if the electrons interact with the same nuclei 
with the same or different strengths. The limit of closely 
spaced electrons is found to support long lived coherence, 
in a manner directly analogous to coherences in photo- 
synthetic complexes (2TJ - Finally, this analysis is used to 
propose a new version of the radical pair model, which 
gives a biologically useful signal in the limit of strong in- 
teractions between the system and the bath, and gives a 
simple explanation of many behavioral experiments. 

The Hamiltonian for a system of electronic spins Si 
and nuclear spins I k in an external field B is given by 
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Here the nuclear magnetic moments are small, so that the 
orientation of nuclear spins due to the weak geomagnetic 
field can be neglected. The orbital angular momentum 
of the radical electrons is assumed to be fixed and to 
contribute no dynamics, while ^2 ■ Si ■ r k = in the m s = 
subspace. Eliminating these terms yields a Hamiltonian 
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where the interaction of each electron with its nearest 
neighbor nuclei (indexed by fcj) has been separated from 
its interaction with the more distant nuclei of the sur- 
rounding bath. 

An electron's interaction VgpjJ • Si with its immediate 
surroundings may rival or exceed the strength of the Zee- 
man interaction B ■ Si. By lifting the degeneracy between 
|tl) and \i-~t), such interactions drive the coherent oscilla- 
tion between singlet and triplet necessary for the radical 
pair model. As the precise receptor or receptors involved 
in the avian compass have not been identified (see [23T - [25j | 
for candidates) , the magnitude of this effect is unknown, 
although the narrow range of frequencies causing disori- 
entation suggests that it is weak (24{. In the absence 
of a known receptor, this paper replaces the Zeeman 
and short range terms with j3 z A/x/2(|t|) (til - lit) (lt|), 
where B z is the magnetic field component in some direc- 
tion preferred by the molecule and Afj, reflects the differ- 
ences in the effective magnetic moment between the two 
electrons once this interaction is taken into account. 

Decoherence between spin states of the radical pair 
occurs due to hyperfine interactions with the surround- 
ing nuclei. Particularly important in this interaction is 
the form of the decoherence - the set of preferred states 
between which decoherence occurs. Both the form and 
the rate of decoherence are determined by the geome- 
try of the excited state. Due to the r -3 dependence of 
the hyperfine interaction, distant electrons will interact 
with effectively distinct reservoirs, while electrons in close 
proximity will have an interaction which may be written 
as the sum of a symmetric term 
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and an antisymmetric term 
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(5) 



where |ri2 is the distance between the electrons, A — 
^{2fi B lf), S = Si + S 2 , AS = 5i - £2 and is the 
angle between ria an d -Rfc- 

If the nucleus is assumed to be in thermal equilibrium 
with a local reservoir, a master equation for the decay of 
electronic coherence due to the hyperfine interaction may 
be found by tracing over the nuclear and environmental 
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degrees of freedom in the density matrix equations of 
motion for the full system. Working in the interaction 
picture, 



Eq. QT] can be evaluated by diagonalizing V. To this 
end, matrix elements of the form / • S can be evaluated 
using the Clebsch-Gordan expansion 



dAt[V T (t), [F'(t-At)X +n+e (t-At)]], \( Jlj2 )JM) = 



d 

w 

(6) 

where Hq +c gives the interaction of the nucleus with the 
bath, Hq is the Zeeman Hamiltonian for the radical pair, 
H = Hg+H^ +e , and H = H + V. Here V can be either 
the symmetric or antisymmetric part of the hyperfine in- 
teraction. The short time evolution of the density matrix 
in Eq. [5] is given by 
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/(t + At) = e 
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Using the Born approximation, the effects of the rad- 
ical pair interacting with the nucleus can be separated 
from the effects of the nucleus interacting with the bath. 
The form of the nucleus's interaction with the local bath 
is unknown, but assumed to drive the system toward a 
Markovian thermal equilibrium, with 
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where {\k}} is some basis preferred by the interaction 
and Pfc the probability of state \k) at thermal equilib- 
rium. Meanwhile, the evolution of the density matrix 
due to H and V can be found using the split opera- 
tor method [26|, so that e tHAt can be approximated by 

e iH At/2 e iVAt e iH At/2 an( j 



p(t - At) = U(At)p(t)U(-At), 



-iff At/2 £ iV At e iH At/2 



(9) 



where [/(At) = e~ 

Substituting Eqs. [5] and |H] into Eq. |H] and tracing 
over environmental and nuclear degrees of freedom will 
now give a master equation for the decoherence due to V. 
For noncommuting Hq and V , the oscillating exponential 
operators of Eq. |6] will result in a time integral contain- 
ing many fourier components, while the rapid decay of 
nuclear spin coherence in Eq. [5] will contribute a delta 
function in time. The time integral can be evaluated by 
noting that 
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so that slowly oscillating terms give a larger contribution 
than rapidly oscillating terms. In the limit of a weak 
magnetic field, e zH ° At « 1 over short times, and Eq. [6] 
becomes 



d 
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and the identity 



\j 1 m 1 j 2 m 2 ) (jim 1 j 2 m 2 \JM) 
(12) 



L\ ■ L 2 — - [J 2 - L\ - L 2 2 \, 



(13) 



where J = L\ + L2 and J 2 ,J Z , L\, and L\ are conserved 
by the interaction. Nonzero matrix elements of / • S are 
given in Table HI 

(l,l;t|/-S|l,l;t) = 1/2 
<l,-l;||/-5|l,-l;|) = l/2 
<l,0;t|/-S|l,l;|) = 1/V2 
(l,l;||/-5|l ! 0;t) = l/v / 2 
<l,0;||/-5|l,-l;t> = 1/V2 
(l,-l;t|/-S f |l,0;|) = l/v / 2 
(1,1;||/-5|1,1;|) = -1/2 
(l,-l;t|/-5|l,-l;t> = -1/2 

TABLE I. Nonzero matrix elements of I • S for states 
\S,ms;mi), from Equations 1121 and 1131 

It can now be seen that the symmetric- and anti- 
symmetric parts of the hyperfine interaction cause de- 
coherence between different sets of preferred states. For 
the symmetric term, S — Si + S2, so that states \s) 
and |t) are eigenstates of the interaction with eigenval- 
ues \s, m) = |0, 0) and |1,0). For the antisymmetric 
AS = Si — §2, states |t4-) and |4-t) are eigenstates with 
eigenvalues \As, Am s ) = |1, 1) and |1, —1). Thus, the ge- 
ometry of the excited state will determine the basis set in 
which off diagonal elements decay most rapidly. Substi- 
tuting V — al ■ S into Eq. [TT] yields an exponential rate 
of decay for the off diagonal elements in the preferred 
basis 



Pi 



' ^ a Pl,mi;2,m,i 



(14) 



Although beyond the scope of the current paper, this 
analysis may be used to describe dynamics in the full 
spin space of the radical pair, as well as the case when 
B ^ 0. In the limit of strong magnetic fields, applicable 
to many laboratory experiments testing the radical pair 
model, the rate of decoherence goes as rather than V, 
while in the zero field limit the m s — subspace may lose 
population due to spin exchange with the surrounding 
nuclei. As ps » pi, the latter process is off resonant 
for even weak magnetic fields and has been neglected. 

The loss of coherence due to the spin bath as a whole 
can be found via a spatial integral assuming a constant 
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density of nuclei per unit volume, so that 

--Rmax 



r = 4tt 



r dr\(r) 



(15) 



and 



Ar = 2tt 



sm(0)dO 



r 2 drAX(r), (16) 



where A(r) = ^Ar~ 3 , AX(r,9) = 2A\r 12 \ cos 6>r" 4 , and 
A=£(2p B ^). 

In Eqs. [15l the formal divergence of the T integral 
is cut off by assuming a slow gaussian decay term in 
addition to the exponential decay from Eq. [14] If 
the decoherence due to a particular nucleus is given by 
/ r (i) = e - A ( r )l*l- <Tt t with a small, then the Fourier 
transform of f r (t) is given by ¥[f r (t)] » F[e" A Wl t l] for 
A(r) >> a and F[/ r (i)] « F[e _<Tt2 ] for a » A(r). Using 
the Fourier convolution theorem, 



F[e 



-T\t\-St 2 



re-MnOltllU 



(17) 
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FIG. 1. (Color online) Decay of coherence in the Bloch sphere 
picture. For initially pure states, trajectories begin on the 
unit sphere before undergoing rapid decay of the z component 
and slow decay of the x and y components. Here, T = 10 
gauss-/is, B — 50 gauss, Afi/fiB = 6 x 10 -4 , Ar = and 
At = 0.1/A+. 



where II is the convolution rather than the multiplica- 
tive product. Noting that F^ 1 [IF[ e - Al 1*1] * F[ e - A2 l*l]] = 
e -(Ai+A 2 )|t| j the effect f the slow gaussian decoherence is 
to introduce a cutoff radius i? max such that A(i? max ) = a, 
beyond which nuclei do not contribute to the exponential 
part of the decoherence. 

For the purpose of modeling the avian compass, it is 
sufficient to model the exponential part of the bath de- 
coherence using the Lindblad equation 



«=| s ),|i),|t4-),|4-t) 



[p], (18) 



where C K [p] = —pL\L K — L^ K L K p + 2L K _pL\ is the Lind- 
blad superoperator corresponding to projection operator 
L K = \k) (k\ and r K = f/2 for \k) = \s) or \t), while 
r K = Ar/2 for |k> = |t|) or ||f). 

The evolution of the density matrix components can 
be found analytically by calculating the evolution due to 
H and T in the singlet /triplet basis, then transforming to 
the |tl), lit) basis to include the effects of Ar. As in |2lj] . 
differential equations for the density matrix components 
due to H and T are given by 



Ptt) 



Ptt) + (B z Ap) 2 ( Pss - Ptt ) =0 



d 2 d , 

J^iPtS - Pst) + 1 foKPtB - Pst) 



(B z Ap) 2 { Pts -p s 



dt 
d_ 
dt 



{pts + Pst) = -r(p ts + p s t) 

{pas + Ptt) = 0. 



Here (p ss — ptt) and {pt s — Pst) behave as damped har- 
monic oscillators, with time dependence 



P(t) = Ae x+t 



Be 



where 



A j 



-f ± ^r 2 - 4(B Z A M )5 



(19) 



(20) 



In the limit that f >> \2B z Ap,\, the system is strongly 
overdamped and the coherence terms will decay much 
more slowly than the base rate of dephasing, T. It is 
this slow loss of coherence which allows for a biologically 
useful signal. 

Transforming to a Bloch sphere representation gives a 
simple picture of the coherence dynamics in this system. 
If \U) -> |0> and ||t) -> 1, then {p ss + p tt ) -> (poo + 
pn) = I, (p ss - ptt) -4 (poi + pio) = xcr x , (p at + pts) -> 
(Pio - Poi) = iyvy, and (p st + p ts ) (poo - Pn) = zcr Zl 
where a x . y . z are Pauli matrices. Decoherence due to the 
Ar term is included by multiplying pio, poi by e~ Ar *. In 
the overdamped limit, (poo — Pn) will decay rapidly as 
e~ rt , so that the z coordinate of a point in the Bloch ball 
will rapidly decay to zero. However, both the x and y 
components are long lived, decaying as e~ x+t . This decay 
is shown grapically in Figure [TJ In spherical coordinates, 
the azimuthal angle 9 is destabilized by interaction with 
the bath while the polar angle <j> is stabilized, so that 
the symmetry group of the long-lived information is U(l) 
rather than SU(2). 

In the radical pair model, creation of a purely singlet 
or purely triplet excited state places the Bloch vector 
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in the |fl), ||t) basis, p = 
overdamped limit, p = —1/2 
1/2 



and in the 



on the equator, where it has no rapidly decaying com- 
ponents. This gives a simple explanation for the long 
coherence times needed to be consistent with disorien- 
tation by a 15 nT oscillating field. However, in con- 
trast to the original radical pair model, <j> does not vary 
with time, so the system does not oscillate between sin- 
glet and triplet character. Rather, it is the decoher- 
ence between states \\\) and 4t) which generates an 
observable signal. For an initially pure singlet state 

' 1/2 1/2 
1/2 1/2 

A+ + Ar 
a+ + Ar o 

-A+ - Ar o 
o a+ + Ar 

sis. Foss of coherence thus manifests itself as a transfer 
of population from the singlet to triplet state, at a rate 
which varies linearly with A + « (B z Ap) 2 /T . 

The dynamics derived above suggest a simple yet ro- 
bust mechanism for a magnetic compass based upon an 
overdamped radical pair. Here the rate of decoherence, 
proportional to \Bq cos(#)| 2 , provides the desired signal 
in the form of population transfer between the singlet and 
triplet states. If the triplet reaction rate is sufficient to 
prevent backwards transfer, the ratio of triplet to singlet 
byproducts is simply 



in the singlet /triplet ba- 



Rts(0) 



A+ + Ar \B Apcos( 



rfc, 



TT- < 21 » 



so that the ratio of triplet to singlet byproducts bears 
a simple relationship to the magnetic field orientation. 
Note that, as this signal is proportional to cos 2 (0), this 
mechanism yields a compass which is sensitive to the ori- 
entation of magnetic field lines but not their polarity, 
consistent with behavioral experiments. 

While the identity of the compass molecule remains 
unknown, some general design considerations may be in- 
ferred. Using Eq. [21] the sensitivity of the compass 
mechanism may be found by calculating the contrast be- 
tween North/South and East/West alignment 



contrast = 



{B Apf 



RtM-R ts {ir/2) = 
R ts (Q) + Rts(ir/2) 2Arr + (BoA/z) 2 



_(22) 

Here the contrast depends upon the ratio of Arr and 
B$Ap. Because the overdamped limit requires that 
r 2 >> (BqAp) 2 , a sensitive compass will require that Ar 
be very small or even zero. Thus, the compass molecule 
may involve an excited state with a geometry which min- 
imizes the antisymmetric part of the hyperfine interac- 
tion, or else a large region where water molecules are 
excluded, so that R m in in Eq. [16] is large. 

A second design consideration relates to the singlet 
reaction rate, k s . As seen in Eqs. I2T1 and l22| k s is a 
free parameter which controls the visibility of the Rt s (0) 
signal without affecting the sensitivity of the compass 



mechanism as a whole. Because the rate of dephasing, 
and triplet formation goes as the square of the magnetic 
field strength, an ability to enhance or inhibit k s would 
allow a migratory animal to use the same biochemical 
pathway for orientation regardless of the local magnetic 
field strength, while a static k s would yield an optimal 
signal at some predetermined field strength. 

The light-activated radical pair compass is a unique 
example of an intrinsically quantum mechanical process 
which survives and even depends upon a strong inter- 
action between the system and the surrounding bath. 
Contrary to expectation, the interplay between coher- 
ent Hamiltonian evolution and bath-induced decoherence 
stabilizes part of the quantum information, allowing an 
organism to sense the orientation of the weak geomag- 
netic field. The avian compass, then, reflects a simple 
model system for the largely unexplored role of quantum 
mechanics in biology. 
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